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Abstract. We formulate a method to compute breakup 
processes of three-body halo systems reacting with a target. 
The reaction between one of the particles and the target as 
well as the corresponding pairwise final state interactions are 
carefully treated. Both absorption and diffraction processes 
are included. Detailed differential and total cross sections are 
calculated for two-neutron (2n) removal from 6 He by fragmen- 
tation on 12 C. Effects due to the core shadowing are taken into 
account. Good agreement with available experimental data is 
obtained. 

25.60.-t - Reactions induced by unstable nuclei 
25.60.Gc - Break-up and momentum distributions 
21.45.+V - Few-body systems 

Introduction. Breakup reactions of three-body systems 
are difficult to describe in general if the formulation is 
required to be both practical and fairly accurate. How- 
ever, due to the steady improvement of radioactive beam 
facilities an increasing amount of fragmentation data de- 
mand analyses where interpretation in terms of few-body 
reactions appear to be an efficient framework jl]-|| . Pre- 
vious attempts are based on the eikonal approximation 
and the adiabatic assumption of frozen intrinsic degrees 
of freedom during the collision [p|-[l3[ . In most of these 
approaches very simple three-body wave functions are 
used and the final state interaction often needs better 
treatment. 

The crucial ingredients in these processes are the fi- 
nal state interaction PJj|,[L4| and the interaction between 
one of the three particles (the participant) and the tar- 
get. The remaining two particles are essentially spec- 
tators and reactions simultaneously involving more than 
one particle are much less important. The large spa- 
tial extension of halo nuclei guarantees that this model 
is a good first approximation which could be improved 
in future more elaborate investigations if needed. We 
shall here fully exploit these simplifying assumptions to 
keep the approach practical enough to compute many ob- 
servables. We shall describe the reaction between each 
particle and the target by the optical model. The final 
state two-body interactions are those used to compute 
the bound state wave function of the three-body projec- 
tile. 

Theory. When a spatially extended three body halo 
hits a relatively small target at high energy the proba- 
bility that more than one of the constituents interacts 
strongly with the target is small. The cross section can 
then be written as a sum of three terms where each 



term describes the contribution to the reaction caused by 
the interaction between the target and the corresponding 
particle 



(i) 



where v is the relative projectile-target velocity, is 

(i) 

the transition matrix element and dv ^ is the density of 
final states. Each of the three terms includes elastic and 
inelastic scattering of the corresponding particle on the 
target. 




FIG. 1. Sketch of the reaction and the coordinates used. The 
target is labelled by and k} label the particles within the 
three-body projectile. 

Neglecting the Coulomb interaction and assuming that 
the target has zero spin we write the transition matrix 
element for the elastic contribution as (compare with the 
formulation for a weakly bound projectile jl5[) 
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where Vch is the interaction between particle i (the par- 
ticipant) and the target, ty JM is the projectile internal 
wave function, and d&^ + ) v are the distorted 

waves in the corresponding two-body subsystems and 
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Xs^, Xs jk x jk are the spin functions (s jk = Sj + s k ). 
The coordinates r.,-fc,R and R' arc dehncd in fig. 1 and 
the conjugate momenta are denoted by the correspond- 
ing p. Primes are used for the final states. The density 
of final states is 



dvf = 5(E' 0i - E 0i ) 
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where is the relative energy of particle i and the 
target and 
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The diffraction cross section for non-polarized projec- 
tile is then obtained after summing over final and aver- 
aging over initial states 
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The absorption cross section where only two projectile 
particles survive in the final state is simply obtained by 
replacing the elastic a e i by a a bs, i-e. 
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The total cross section arising from Voi is the sum of 
eqs.|]) and (pi]). The function M^.^ s . is the transi- 
tion matrix computed in ]ll| . 

Parameters. We now consider the nucleus 6 He 
(n+n+a) with the wave function obtained by solving the 
Faddeev equations in coordinate space [jl0| by using the 
potentials from [ fl6[ . The resulting three-body wave func- 
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tion has 88% of jr- and 12% of s 2 -configurations. The 
binding energy is 0.95 MeV and the root mean square 
radius is 2.45 fm. 

For the neutron-target interactions we use non- 
relativistic optical potentials obtained from relativistic 
potentials through a reduction of the Dirac equation into 
a Schrodinger-like equation jlTJ]. In particular we focus 



on a carbon target, and for the neutron- 12 C interaction 
we use the parametrization EDAI-C12 [H valid for a 
range of neutron energies from 29 to 1040 MeV. We in- 
clude 35 partial waves in the calculations. 

The detected 2n removal contribution from the inter- 
action between a-particle and target is expected to be 
relatively small. First the elastic scattering cross section 
for a-particles on 12 C is about 35% of the total cross 
section |19[ ], which in turn is of the same order as the 
neutron C total cross section. Second a substantial 
fraction of these 35%, where the a-particle survives the 
reaction with the same energy, are in fact elastically scat- 
tered 6 He-particles and as such not contributing to the 
2n removal cross section. Third the very forward angles, 
containing most of the elastic cross section, are experi- 
mentally excluded. We shall therefore neglect this con- 
tribution in the present letter. 

The final state two-body continuum wave functions 
(<ft> k ( + \ (f> 0l ( + ') are now calculated with the appropri- 
ate boundary conditions and normalization |11|. Then 
eqs.(||) and ( |ll|) are integrated over the unobserved quan- 
tities and different momentum distributions or differen- 
tial cross sections are obtained. However, the so-called 
core shadowing problem remains p0|~p2]|. Basically this 
means that the finite sizes of the target and the a-core 
prohibit some of the neutron-target reactions namely 
those where the a-particle also interacts with the target. 
In the eikonal approximation this effect is accounted for 
by multiplying the halo wave function by a profile func- 
tion. A simpler approximation is the black disk model 
where the profile functions are step functions pH , and 
the interior part of the halo wave function is eliminated. 
This amounts to eliminate the configuration which is not 
consistent with the reaction in question. 

In particular, in a process where the neutron is re- 
moved by the target without destroying the a-core, we 
must exclude the part of the three-body wave function, 
where the removed neutron is too close (or even inside) 
the core. We do this by assuming that 'F; = when the 
distance between the interacting neutron and the center- 
of-mass of the remaining two projectile particles is suf- 
ficiently small. This amounts to assuming that the Ja- 
cobi coordinate y = r n ^ an y/ (m a + m n )/ (m a + 2m n ) is 
smaller than a certain cut off value y c , which presumably 
is better left as a parameter, but clearly is related to the 
sizes of the reacting particles. In fact, the sum of the 
target and core radii is around 4 fm, and we expect that 
for r„ !Cm < 4 fm the neutron can not be removed without 
disturbing (or breaking) the core. 

Observables. In fig. 2 we compare the calculated and 
experimental [^) transverse 5 He momentum distribution 
after 6 He fragmentation on 12 C at 240 MeV/u. Without 
shadowing a too broad momentum distribution is found, 
while the shadowing with y c = 6 fm gives a too narrow 
distribution. For the estimated value of y c — 4 fm, the 
central part of the distribution up to about 60 MeV/c 
is rather well reproduced. The momentum tails receive 
contributions from distances where the neutron is inside 
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the core. The model is in other words not applicable for 
these momenta. Thus y c = 4 fin seems to be an appro- 

friate value. This is consistent with the interpretation in 



1000 



CL 



800 



600 



400 - 



200 




-100 



p t ( He, MeV/c) 



FIG. 2. Transverse 5 He momentum distribution after fragmen- 
tation of 6 He on 12 C at 240 MeV/u. 

The absolute value of the 2n removal cross section is 
sensitive to the value of y c as seen in fig. 3a, where 
we show results for breakup of 6 He on 12 C as function of 
beam energy for different shadowing parameters. The ab- 
sorption processes amounts to 70-75% of the cross section 
depending on energy, consistent with the known data of 
neutron scattering on 12 C iQ. In all cases the shadowing 
for y c — 4 fm reduces the cross section by a factor around 
3.6, but the absolute cross sections should be taken with 
caution due to the high sensitivity to y c . However, the 
numbers are encouraging and in agreement with both 
experimental values 23 and theoretical estimates PJl^f. 
The predicted minimum at about 200 MeV arises from a 
corresponding minimum in the n- 12 C cross section [ Q . 

In fig. 3b we show the value of the total 2n removal 
cross section, with y c equal to 3, 4 and 5 fm relative to its 
value at a beam energy of 100 MeV/u (0.430 b for y c = 3 
fm, 0.250 b for y c = 4 fm and 0.130 b for y c = 5 fm ). 
We also show the width of the transverse neutron mo- 
mentum distribution relative to the width at 100 MeV/u 
(77 MeV/c, 75 MeV/c and 73 MeV/c for y c =3, 4 and 5 
fm, respectively). Both these relative quantities are es- 
sentially independent of the shadowing. In addition we 
note that the width of the neutron momentum distribu- 
tion is almost constant for these high beam energies. 

In fig. 4a we compare the calculated and measured 
transverse neutron momentum distribution after 6 He 
fragmentation on carbon. 
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FIG. 3. (a) Calculated 2n removal cross section after 6 He frag- 
mentation on 12 C as a function of the beam energy, (b) The same 
reaction as (a). The circles and squares respectively indicate the 
ratio between the total 2n removal cross section and the width of 
the transverse neutron momentum distribution as a function of the 
beam energy. 

We also show the results of calculations with and with- 
out final state interaction based on the sudden approx- 
imation as described in 10 nj . These results resemble 
those without shadowing. This fact supports the validity 
of the sudden approximation with inclusion of final state 
interaction as a first approximation to the description of 
these reactions. However, even if this approximation is 
rather good, the agreement with the experiment is im- 
proved when shadowing with y c = 4 fm is used. 
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FIG. 4. (a) The transverse neutron momentum distribution af- 
ter 6 He fragmentation on 12 C at 240 MeV/u. The results of the 
calculations described in ref. [11] with (long dashed curve) and with- 
out (short dashed curve) inclusion of final state interactions, (b) 
Computed longitudinal and transverse neutron momentum distri- 
bution for the same case as (a). 
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In fig. 4b we compare the calculated longitudinal and 
transverse neutron momentum distributions. They are 
very similar, but a small asymmetry, due to the energy 
dependence of the n- 12 C optical potential, is present in 
the longitudinal distribution. 

Conclusions. A practical method treating one particle- 
target interaction at a time is formulated for spatially 
extended three-body systems colliding with a light tar- 
get. Absolute values of differential cross sections includ- 
ing their dependence on beam energy and target struc- 
ture can then be computed. The longitudinal and trans- 
verse momentum distributions in breakup reactions are 
now distinguishable. Both diffractive and absorptive pro- 
cesses are treated.. The method is subsequently applied 
to breakup reactions of the halo nucleus 6 He (n+n+a) on 
a 12 C target. The final state interaction and the shadow- 
ing are crucial for the neutron momentum distributions 
and the absolute cross sections respectively. Remark- 
able agreement is achieved with the observed quantities 
for absolute as well as for differential cross sections. If 
not a correct description we provide at least a very good 
parametrization. 
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